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Abstract 

Extraterrestrial neutrinos can initiate deeply developing air showers, and those that traverse 
the atmosphere unscathed may produce cascades in the ice or water. Up to now, no such 
events have been observed. This can be translated into upper limits on the diffuse neutrino 
flux. On the other hand, the observation of cosmic rays with primary energies > 10 10 GeV 
suggests that there is a guaranteed flux of cosmogenic neutrinos, arising from the decay of 
charged pions (and their muon daughters) produced in proton interactions with the cos- 
mic microwave background. In this work, armed with these cosmogenic neutrinos and the 
increased exposure of neutrino telescopes we bring up-to-date model-independent upper 
bounds on the neutrino-nucleon inelastic cross section. Uncertainties in the cosmogenic neu- 
trino flux are discussed and taken into account in our analysis. The prospects for improving 
these bounds with the Pierre Auger Observatory are also estimated. The unprecedented 
statistics to be collected by this experiment in 6 yr of operation will probe the neutrino- 
nucleon inelastic cross section at the level of Standard Model predictions. 



*On leave from Institute for Theoretical Physics, Eotvos University, Budapest, Hungary. 



1 Introduction 



Ultrahigh energy cosmic neutrinos may reveal aspects of nature hidden to us so far. They can 
point back to very distant sources, resolving the origin of the highest energy cosmic rays and the 
underlying acceleration mechanism. Upon their arrival at Earth, they interact with nucleons at 
centre-of-mass energies around several hundreds of TeV, probing the energy regime far beyond 
the reach of terrestrial colliders. 

The detection of ultrahigh energy cosmic neutrinos is a very challenging task due to their weak 
interactions. Although the neutrino-nucleon inelastic cross section increases with energy, the cos- 
mic neutrino flux falls even more steeply with energy. Large-scale projects and novel techniques 
are being deployed and proposed. The Pierre Auger Observatory (PAO) [1], followed by Ice- 
Cube [2], and possibly EUSO [3] and OWL [4], will reach a sensitivity of the level of theoretical 
predictions for the cosmic neutrino fluxes. Moreover, before the next generation experiments lead 
us into the exciting era, the non-observation of neutrino-induced events reported by the Fly's 
Eye [5], the AGASA [6,7] and the RICE [8,9] collaborations can be used to set upper limits on 
the diffuse neutrino flux. Additionally, by exploiting a certain prediction for the neutrino flux, 
the search results can also be turned into upper bounds on the neutrino-nucleon inelastic cross 
section [10-12]. 

Ultrahigh energy cosmic neutrinos from diverse sources are predicted, and their existence is 
strongly supported by the observation of the ultrahigh energy cosmic rays (UHECRs). Among 
all, the so-called cosmogenic neutrinos [13] are almost guaranteed to exist. They originate from 
the decay of charged pions produced in the interactions of protons with the cosmic microwave 
background (CMB). There are still some uncertainties in the estimate of the cosmogenic neutrino 
flux. This is due to our poor knowledge of the nature and the origin of the UHECRs. Possible 
ranges for the size of the cosmogenic neutrino flux have been investigated elsewhere [14] assuming 
that the cosmic ray spectrum beyond 10 8 GeV is dominated by extragalactic protons with an 
isotropic distribution of sources. In this work, armed with these cosmogenic neutrinos and the 
increased exposure of neutrino detectors, we show that the neutrino-nucleon inelastic cross section 
is tighter constrained than previously noted. Uncertainties in the cosmogenic neutrino flux are 
discussed and taken into account in our analysis. The prospects for improving these bounds with 
the PAO are also estimated. 

The structure of the paper is the following. In Sec. 121 we examine acceptances for neutrino 
detection. In particular, we compute the effective apertures of AGASA and RICE as examples of 
ground arrays and under-ice neutrino telescopes, respectively. We also estimate the aperture of 
the PAO in order to investigate future sensitivity to physics beyond the Standard Model (SM). 
In Sec. H3 we present an overview of cosmogenic neutrino fluxes. In Sec. 0] we derive model- 
independent upper bounds on the neutrino-nucleon inelastic cross section from the search results 
reported by the AGASA and the RICE collaborations. In Sec. the prospects for improving 
these bounds with the PAO are estimated. Section El comprises our conclusions. 

2 Search for Ultrahigh Energy Cosmic Neutrinos 

Due to their feeble interaction and flux, ultrahigh energy cosmic neutrinos are extremely difficult 
to detect. To overcome this difficulty one either needs a huge target volume for the detectors, 
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or one has to employ novel detection techniques. Besides, the unwanted strong backgrounds 
from cosmic ray protons must be effectively reduced. Following these considerations, ultrahigh 
energy cosmic neutrinos are searched for in the Earth atmosphere [5-7], in the Greenland [15] 
and Antarctic [16] ice sheet, in the sea/lake [17], or even in the regolith of the moon [18]. For 
this purpose fluorescence detectors, ground arrays, underwater and -ice neutrino telescopes exist 
in various stages of maturity, from proposed to nearly completed. In this section we summarise 
general formulae for estimating neutrino-induced event rates. We discuss separately their applica- 
bility in the case when neutrinos interact significantly more strongly than in the SM, as predicted 
in many new physics scenarios [19]. 

We start with the differential rate of air showers initiated at the point (£, 9) by particles (cosmic 
ray protons, cosmic neutrinos etc.) incoming with energy E and of flux F(E), where I is the 
distance of this point to the detector measured along the shower axis, and 9 is the angle to the 
zenith at the point where the shower axis hits the Earth's surface (cf. Fig. The number of 
air showers induced due to the interaction, the cross section for which is cr(E), per unit of time 
t, area A, solid angle Q (with dfl = sin 9 d9 d<p) and energy E^ is 

d / d 4 N \ 1 <7 tot (E)x(e,o) 

where m p is the proton mass, and p a j r is the air density at the altitude h. The relation of the energy 
deposited in a visible shower, E sh , to the incident particle energy E is dependent on the scattering 
process. For protons, E B ^ « E, whereas in general, E^ = yE, where y is the inelasticity. For 
neutrino interactions, the inelasticity distribution is measured through deep inelastic scattering 
processes for which y = [E v — E[)/E Vl where E[ is the energy of the final state lepton. For 
simplicity, throughout this paper we take the average value of this distribution, denoted by (y). 
The exponential term in Eq. ((TJ accounts for the flux attenuation in the atmosphere, where the 
total inelastic cross section a tot (E) = a (E) + a new (E) can receive contributions from SM and 
new physics interactions. 1 

After carrying out the integral of p a i r {£, 9) d£ = — dx(£,9) over the range of depths, X Q b 8 (9) = 
X{9) — X uno (9), within which showers induced are visible to the ground array detectors, the rate 
of neutrino-induced events at a ground array with threshold energy E t ^ can be well estimated 
as [12] 

d6 S(E sh , 6) 2vr sin 6 le 
~ tA p I dE sh ^ ^ F V {E V ) P(E Bh ) att(^) , (2) 

JE th VvNK&v) 

where X(9) is the atmospheric slant depth of the ground array, X uno (6) is the minimum atmo- 
spheric depth a neutrino must reach in order to induce an observable shower to the ground array 
(cf. Fig.HJ), and 

&tt(E„) = / d cos 9 2tt cos 9 e ■ I 1 — e m p ) . (3) 

J cos jW V / 



1 In general, the inelasticities of the Standard Model and of the new physics contributions might be different. 
For cr ncw 3> cr SM , as considered in our analysis, one is sensitive only to the (y) of the new physics. 
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We approximated the effective area as S(E, 9) ~ A p P(E) cos 9, with A p a parameter which is of 
the same order of the detector's geometric area. 

In the SM, the hadronic component produced in the charged current (i/N —>■ IX) and neutral 
current neutrino-nucleon interactions (i/N — > uX) generates a shower which is visible at the 
ground array. The hadronic component inherits an energy E s ^ = (y) E v from the neutrino. The 
charged current process for electron neutrinos produces both a hadronic shower (E s h ~ 0.2E U ) 
and an electromagnetic shower (E S h « 0.8E U ). Thus, in such a case all the energy is deposit in 
the atmosphere. 

Note that the neutrino-nucleon inelastic cross section a V N{E u ) is also contained in the (over zenith 
angle integrated) "attenuation factor" a,tt{E v ). This factor determines the effective aperture for 
neutrino detection. In the SM, neutrinos interact weakly (e.g. a^Jf ~ 10~ 4 mb at E v = 10 11 GeV), 
so Eq. (J2J) can be reduced to the more familiar formula 



It is obvious that ground arrays should look at the quasi-horizontal direction, i.e. 9 > 70°, in 
order to achieve a large range of observability X Q b s , and thus a larger neutrino detection rate. At 
the same time, it helps to reduce the backgrounds from cosmic ray protons. 

For an under-ice neutrino telescope, the hadronic and electromagnetic cascades are detectable 
only when they are initiated within or very close to the detector. The rate for the contained 
events per solid angle Q is 



where p lce = 0.92 g/cm 3 is the ice density, and V is the effective volume of the detector, which 
is usually energy dependent. The exponential function takes into account the neutrino flux 
attenuation in the atmosphere and in the ice. Here, x{9) = X atm (9) + p icc d(9), where d is the 
distance traversed by neutrinos incident with zenith angle 9 (measured from the centre of the 
detector) from the ice upper surface to the detector's surface. 

The flux attenuation in the ice is much more effective than in the atmosphere. The number of 
neutrinos reaching the detector's target volume will be strongly reduced if the neutrino-nucleon 
cross section is enhanced. On the other hand, at sufficiently high energies, under-ice neutrino 
telescopes are also sensitive to interaction vertices located several kilometers away which produces 
muon(s) or tau(s) passing through the detector (through-going muon and tau events). This 
increases their effective aperture for neutrinos largely. However, beyond the SM, the predictions 
of the muon or tau spectra depend strongly on the scenario considered (see e.g. [20]). 

Tau neutrinos can also be searched for through the showers induced by tau lepton decay in the 
atmosphere [21,22]. Though conventional astrophysical sources do not produce tau neutrinos, 
terrestrial experiments (see e.g. [23]) have shown that and v T are maximally mixed with a 
mass-squared difference ~ 10 -3 eV 2 . This, together with the known smallness of | 1 2^3) | 2 [24], 
implies that the is^'s will partition themselves equally between z/^'s and ^ T 's on lengths large 
compared to the oscillation length A osc ~ 1.5 x 10~ 3 (E v /PeV) pc; here u 3 = {y^ + i/ T )/y2 is 
the third neutrino eigenstate. One consequence of this remarkable symmetry is the process of 





(5) 
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Figure 1: Schematic diagram for our definition for the range of depths X ^ s (8) for a ground array 
situated at an atmospheric slant depth X{6). Showers initiated below X mio (9) are visible to the 
ground array. On average, a 10 10 GeV shower develops to its maximum X max after traversing 
about 800 g/cm 2 in the atmosphere. 



decohereing: any initial flavour ratio w a (a = e, /i, r) that contains w e = 1/3 will arrive at Earth 
with equipartition on the three flavours [25]. Since cosmic neutrinos originate dominantly from 
the decay of 7r =1= and their muon daughters, their initial flavour ratios w e : : w T = 1 : 2 : 0, 
should arrive at Earth democratically distributed, i.e., 1:1:1. 

If the geometrical configuration is met, tau neutrinos entering the Earth at large nadir angles will 
convert to a tau lepton near Earth's surface. The tau lepton exits the Earth and then decays in 
the surroundings of the array. The detection rate depends on the neutrino survival probability 
after a distance z in the Earth for conversion to a charged lepton, and on the probability the 
charged lepton would escape the Earth's surface with an energy above the detector's threshold, 
i.e. [22] 

oc / dz e- J ° dz " N[ u)Ps P esc (E T ) ■ (AQU(E T ) , (6) 



dt dE T J m p 

where p(r) is Earth density at radius r, and p s ~ 2.65 g/cm 3 Earth's surface density. The effective 
aperture (AQ) e s depends on the probability of the r inducing a shower which triggers the detector. 

While improving the prospects of neutrino detection to a large extent, this signal is not always 
present when neutrino interactions predicted by scenarios beyond the SM are considered. An 
enhancement of the neutrino-nucleon total cross section will lead to a decrease of the neutrino 
interaction length in the Earth. This will make the geometrical configuration even more difficult to 
meet. Though the neutrino interaction probability is enhanced, the hadronic interaction products 
will be absorbed in the Earth. The leptons which can escape the Earth may only inherit a tiny 
fraction of the incident neutrino energy, so that the shower they induce is practically undetectable. 
This has been demonstrated in [26] with the example of black hole production in ultrahigh energy 
neutrino-nucleon scattering. Nevertheless, the earth-skimming tau neutrinos provide a way to 
discriminate new physics interactions from SM "backgrounds" . 

Up to now, experiments with ultrahigh energy neutrino sensitivity include AGASA [6,7], RICE [8, 
9], FORTE [15], GLUE [18], Fly's Eye [5], and AMANDA [16]. In the following we consider 
AGASA and RICE as case studies: i) AGASA is sensitive to deeply developing showers in the 
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atmosphere and uses a well developed technique; ii) RICE has a greater exposure, though the 
experimental technique has a shorter track record. In the future, the best measurements will 
come from the PAO [27] and IceCube [28] experiments, which have a similar aperture at these 
energies; below we consider the PAO as an example of future sensitivity to physics beyond the 
SM. 

2.1 AGASA 

The Akeno Giant Air Shower Array (AGASA) occupies farm land near the village of Akeno 
(Japan) at a longitude of 138°30' East and latitude 35°30' North [29]. The array, which consists 
of 111 surface detectors deployed over an area of about 100 km 2 , has been running since 1990. 
About 95% of the surface detectors were operational from March to December 1991, and the 
array has been fully operational since then. A prototype detector operated from 1984 to 1990 
and has been part of AGASA since 1990 [30]. 

The AGASA Collaboration has searched for deeply-penetrating inclined (zenith angle 6 > 60°) 
air showers. Non-observation of neutrino events during a running time of 9.7 x 10 7 s was reported 
in [6], implying, for zero events of background, an upper limit of 2.44 at 90% confidence limit 
(CL) [31]. The AGASA data for deeply penetrating air showers recorded from December 1995 to 
November 2000 (with an effective lifetime of 1710.5 days) is published in [7]. Deeply penetrating 
events must satisfy following search criteria: i) 9 > 60°, ii) |A^ ax — X s i ant | < 500 g/cm 2 , where 
Aslant is the atmospheric slant depth of the AGASA array centre along the shower axis, and Hi) 
A^ ax , A^ ax > 2500 g/cm 2 , where r\ and 5 parametrise the lateral distribution of shower particle 
densities and the curvature of the shower disc front, respectively. By fitting them to the empirical 
formulae, the shower maximum A max can possibly be deduced. Specifically, there was one event 
observed, consistent with the expected background from hadronic showers l^t^l^ (MC 
statistics and systematic). The AGASA search result therefore corresponds to an upper bound 
of 3.5 events at 95% CL [31]. In Ref. [11], the AGASA search result has been combined with the 
non-observation of deeply-penetrating showers reported by the Fly's Eye Collaboration [5] to set 
model independent upper limits on the ultrahigh energy cosmic neutrino flux (cf. Fig. |2J). 

The effective aperture for deeply penetrating showers has been parametrised in [12]. In short, 
we set X uno (9) = X{6) — 1300 g/cm 2 and X uno > 1700 g/cm 2 , in order to take into account 
the fact that on average a 10 10 GeV shower traverses about 800 g/cm 2 before it develops to its 
maximum [32]. The detection efficiency P(E s h) reaches 100% at w 10 10 GeV, where the 
aperture for electromagnetic showers induced by u e is found to be A e s ~ 300 m 2 sr [7]. The event 
rate at AGASA is given by 

^ = J dE sh F(E u )A cS (E u ) , (7) 

where 

Aes(E v ) = A p P(E sh ) att(£„) . (8) 

From Eq. (JHJ) one can reliably deduce the effective area A p . For u e electromagnetic showers we 
obtain Ap™ = 56 km 2 , with E^ = E v . As there is no estimate of the aperture for hadronic showers 
available, to be conservative hereafter we simply assume A^ A = A™. At energies < 10 10 GeV, 
we interpolate between P(E sh ) = at 10 s GeV and P(E sh ) = 1 at 10 10 GeV. 
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2.2 RICE 



The Radio Ice Cherenkov Experiment (RICE) at the South Pole aims to detect electron neutrinos 
of energies > 10 6 GeV based on the principle of "radio coherence" . Namely, the electrons produced 
in the v e charged current interactions {v e N — > e~X) in the Antarctic ice-cap initiate sub-showers 
which emit Cherenkov radiation over a wide range of electromagnetic frequencies. The RICE 
detector is sensitive only to the long-wavelength range. Monte Carlo simulations show that RICE 
achieves an effective volume V e s ~ 1 km 3 at E s ^ m 10 8 ' 5 GeV and approaches V e g ~ 10 km 3 
at higher energies [8]. Furthermore, since hadronic showers do not suffer significantly from the 
Landau- Pomeranchuck-Migdal (LPM) effect [33], they can increase the detection efficiency of 
an under-ice radio Cherenkov detector at energies in excess of 10 6 GeV. As no candidates for 
neutrino-induced events have been seen during data-taking in 1999, 2000 and 2001 (3500 hrs 
lifetime), 95% CL upper limits on the diffuse neutrino flux have been derived [9]. 

In order to calculate the expected event rates with Eq. we approximate the effective target 
volume of the RICE detector as follows. Due to its location in ice (at 100 - 300 m depths), there 
are no hadronic backgrounds, so RICE utilises the whole zenith angular range 9 = 0° — 90°. In our 
approach, the Monte Carlo effective volume for the electromagnetic and the hadronic cascades 
given by the RICE Collaboration in [9] is approximated as V e g(E) « irr 2 (E) z(E). The height of 
this cylinder is fixed to z — 1 km from the ice surface to below, since the most detection efficiency 
is in the upper km of the ice [8] . The zenith angle 9 is measured from the centre of the detector, 
i.e. at a depth of 0.5 km. With the knowledge of the detector geometry we are able to estimate 
the attenuation of the neutrino flux in the ice properly. 

We have checked that, replacing our approximation for the RICE effective aperture in Eq. (JSJ), 
we were able to reproduce the RICE upper limits on v e flux reported in [8] to within « 35%. 
This may be further reduced (to ~ 20%) if one takes into account nuclear structure effects on the 
neutrino- nucleon charged current cross section [34], as was done in Ref. [8]. 

Following the procedure introduced in Ref. [11], one can easily update the model- independent 
upper bounds on neutrino fluxes using the search result reported by the RICE Collaboration [9]. 
These new limits, which are shown in Fig. |2J surpass the previous estimate derived by combining 
Fly's Eye + AGASA exposures by about one order of magnitude. 

2.3 PAO 

The Pierre Auger Observatory [1], which is actually comprised of two sub-observatories, will be 
the next large-scale neutrino detector in operation. The Southern site is currently operational 
and in the process of growing to its final size of A ~ 3000 km 2 . Another site is planned for 
the Northern hemisphere. The PAO works in a hybrid mode, and when complete, each site will 
contain 24 fluorescence detectors overlooking a ground array of 1600 water Cherenkov detectors. 
During clear, dark nights, events are simultaneously observed by fluorescence light and particle 
detectors, allowing powerful reconstruction and cross-calibration techniques. The first analyses of 
data from the PAO are currently underway [35] and it is expected that first results will be made 
public in the Summer of 2005 at the 29'th International Cosmic Ray Conference. 

For standard neutrino interactions in the atmosphere, each site of PAO reaches about 1.3 x 10 7 kT 
sr of target mass around 10 10 GeV [27]. The sensitivity of PAO for neutrino-induced (i.e. deeply 
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penetrating) hadronic showers, defined as one event per year and per energy decade, is shown in 
Fig. |21 (dashed-dotted line). An even greater acceptance [36] should be achievable for the case of 
Earth-skimming neutrinos which produce a r, that decays and generates a shower in the ground 
array. The projected sensitivity for u T is also shown in Fig. |21 (hatched area). The prospect of 
tau neutrino detection by the PAO fluorescence detector has been also investigated [37], though 
of course the acceptance is decreased because of the 10% duty cycle. 

The rate of neutrino-induced events at the ground arrays of PAO can be calculated using Eq. (J2J). 
We estimate the effective aperture for neutrinos, 

Aff(^) = att (^) p (^) ME sh ) , (9) 

through a comparison with the geometric acceptance given in Ref. [27], where A e s{E u ) »s <j v n(E v ) x 
acceptance (E s -^)/m p . The detection efficiency is found to be 

P(E sh ) = 0.654 log 10 [(£ sh /l GeV) x 10 9 ] - 10.9 , (10) 

for E^ < 10 8 ' 9 GeV, and P(i2 s h) = 1 above this energy. The parameter A p in Eq. (jHJ) is energy- 
dependent, because the PAO acceptance does not saturate in the entire energy range. Our 
selection criteria are i) 75° < 9 < 90° for the zenith angle, and ii) X max > 2500 g/cm 2 for 
the shower maximum, which corresponds to requiring X uno > 1700 g/cm 2 in our approach. We 
consider showers with axis falling in the array. The altitude of the PAO southern site (1200 m 
above sea level) is also taken into account. We obtained Ap ad (E) rs 1.475 km 2 (E/l eV) ' 151 
for hadronic showers above » 10 10 GeV, and A™(E) « 7.037 x 10 6 km 2 (E/l eV) -0 ' 208 for 
electromagnetic showers. As shown in Ref. [27], the aperture for all showers (i.e. including 
showers with axis not going through the array) is roughly 1.8 to 2.5 times larger. 



3 Ultrahigh Energy Cosmic Neutrino Fluxes 

The opacity of the CMB to ultrahigh energy protons propagating over cosmological distances 
guarantees a cosmogenic flux of neutrinos, originated via the decay of charged pions produced in 
the proton-photon interactions [13]. The intermediate state of the reaction p + 7cmb —>■ N + tc 
is dominated by the A + resonance, because the n decay length is smaller than the nucleon mean 
free path on the relic photons. Hence, there is roughly an equal number of ir + and tt°. Gamma 
rays, produced via tt° decay, subsequently cascade electromagnetically on the cosmic radiation 
fields through e + e~ pair production followed by inverse Compton scattering. The net result is a 
pile up of 7 rays at GeV energies, just below the threshold for further pair production. On the 
other hand, each 7r + decays to 3 neutrinos and a positron. The e + readily loses its energy through 
synchrotron radiation in the cosmic magnetic fields. The neutrinos carry away about 3/4 of the 
7r + energy, and therefore the energy in cosmogenic neutrinos is about 3/4 of the one produced in 
7-rays. 

The normalisation of the neutrino flux depends critically on the cosmological evolution of the 
cosmic ray sources and on their proton injection spectra [38]. Of course, the neutrino intensity 
also depends on the homogeneity of sources: for example, semi-local objects, such as the Virgo 
cluster [39], could contribute to the high energy tail of the neutrino spectrum. Another source 
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of uncertainty in the cosmogenic neutrino flux is the Galactic — > extragalactic crossover energy 
of cosmic rays: while Fly's Eye data [40] seem to favour a transition at 10 10 GeV, a recent 
analysis of the HiRes data [41] points to a lower value ~ 10 9 GeV. This translates into different 
proton luminosities at sources and consequently different predictions for the expected flux of 
neutrinos [42]. 

Very recently, some of us (FKRT) have performed an investigation of the actual size of the 
cosmogenic neutrino flux [14]. The assumptions made therein were i) all observed cosmic ray 
events are due to protons, and ii) their sources are isotropically distributed in the universe. It was 
further assumed that all sources have identical injection spectra J p , and that the redshift evolution 
of the source luminosity and of the source co-moving number density can be parametrised by a 
simple power-law. Therefore, the co-moving emissivity of protons injected with energy Ei at a 
distance r from Earth can be written as 

C p (r, = p (1 + z(r)) n Q(z - z min ) 6(> max - z) J p {Ei) , (11) 

where the redshift z and the distance r are related by dz = (1 + z) H{z) dr. The Hubble expansion 
rate at a redshift z is related to the present one H through H 2 (z) = H^ [fi M (l + zf + tt A ], where 
= 0.3 and Q\ = 0.7 were chosen. The results turned out to be rather insensitive to the precise 
values of the cosmological parameters within their uncertainties. The redshift evolution index n 
accounts for the evolution of source emissivities in addition to the pure redshifting (n = 0). The 
minimal and maximal redshift parameters z m \ n and z max exclude the existence of nearby and 
early time sources. The values z m [ n = 0.012 (corresponding to r min = 50 Mpc) and z max = 2 were 
chosen. 

The propagation of protons towards Earth can be well described [43] by the propagation function 
Pb\ a (E; Ei,r) (here b = a = p). It specifies the probability of detecting a particle of species b 
above an energy E on Earth due to one particle of species a created at a distance r with energy 
Ei. To simulate the photohadronic processes of protons with the CMB photons, the SOPHIA 
Monte-Carlo program [44] was adopted. For e + e~ pair production, the continuous energy loss 
approximation was used. With the help of the propagation function, the number of cosmic ray 
protons or cosmogenic neutrinos (b = v) arriving at Earth with energy E per units of energy, 
area, time and solid angle can then be calculated by 



rp ( rp\ d*N b 1 r° [°° ^ dP blp (E;Ej,r) 

Fb{E) = dEdAdtdn = ^J dEi J drH dE c ^ Ei) - (12) 



This formula can also be easily generalised to arbitrary source emissivity. 

To determine the "most probable" cosmogenic neutrino flux F V (E) from the observed cosmic ray 
spectrum, an E~ a power-like injection spectrum for the protons was assumed in Eq. (jllj) up to 
E max , the maximal energy attainable through astrophysical acceleration processes in a bottom-up 
scenario, i.e., 

J p {Ei) = J E7 a 6(£ max - Ei) . (13) 

The predicted differential proton flux at Earth F P (E) was compared with the most recent observa- 
tions by the AGASA [45] and the HiRes [46] experiments separately. A fitting procedure yielded 
the most probable values for the maximal proton injection energy E may[ , the power-law index a, 
and the redshift evolution index n. The factors Jo and po from Eq. ()11[) served to normalise the 
predicted proton flux. 
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For each E may[ , the compatibility of various (a,n) pairs with the cosmic ray data in the energy 
range between E_ = 10 8,2 GeV and E + = 10 11 GeV was checked. This procedure gave the 
2a confidence regions in the a — n plane for each -E max . The best fit values are found to be 
£ max = 10 12 - 5 GeV, a = 2.57, n = 3.30 for AGASA, and £ max = 10 125 GeV, a = 2.50, n = 3.80, 
for HiRes. The corresponding cosmogenic neutrino flux from fitting to AGASA data is shown in 
Fig. |21 (solid line). The one from fitting to the HiRes data is smaller by roughly a factor 1.1 to 
1.25, with larger 2a uncertainties. 

It should be noted that when performing the compatibility check of (a,n) with E_ = 10 9,5 GeV 
(and E + = 10 11 GeV) for each E max , the redshift evolution index n is no longer constrained 
in this case. To consider the possibility of a Galactic — > extragalactic transition in agreement 
with the Fly's Eye data, in our analysis we adopt the cosmogenic neutrino flux estimates of 
Protheroe and Johnson (PJ) [48]. This analysis incorporates the source cosmological evolution 
from estimates [49] of the power per comoving volume injected in protons by powerful radio 
galaxies. Here we use PJ + estimate with an injection spectrum with -E max = 10 12 ' 5 GeV. 
We stress that the PJ flux agrees with a most recent estimate [50] in the entire energy range, 
whereas the spectrum obtained in earlier calculations [38] is somewhat narrower, probably as a 
result of different assumptions regarding the propagation of protons. 

Cosmogenic antineutrinos can also be produced via decay of neutrons photo-dissociated from 
heavy nucleus primaries by the CMB. However, it turns out that antineutrinos from neutron 
/5-decay contribute relatively little to the cosmogenic flux in the energy region of interest (see 
Appendix and Fig. |2J). 

4 Upper Bounds on the Neutrino-nucleon Cross Section 

Upon their arrival at Earth, ultrahigh energy cosmic neutrinos interact with nucleons at centre- 
of-mass energies \fs approaching several hundreds of TeV. They probe the energy regime far 
beyond the reach of terrestrial colliders. In the SM, ultrahigh energy neutrinos scatter deep- 
inelastically on nucleons. The double differential DIS-cross section d 2 a/dxdQ 2 can be calculated 
with the help of the structure functions, where x and — Q 2 are the Bjorken variable and the 
invariant momentum transfer, respectively. The kinematic region in (Q 2 ,x) probed thereby is 
Q 2 ~ m 2 w ~ 6.4 x 10 3 GeV 2 and x = Q 2 /(ys) w 1.7 x 10' 7 /(E u /10 u GeV), not accessible by 
HERA or other collider experiments. It is therefore necessary to reliably extrapolate the nucleon 
structure functions to the region of high Q 2 and very small x values. The cross sections for the 
uN charged and neutral current interactions have been estimated by several approaches [51-54]. 
All predict a power-like growth behaviour with energy (cf. Fig. |3]for one of the predictions [51]). 
This reflects simply the rapid increase of parton densities towoards small Bjorken x, as predicted 
by perturbative QCD and confirmed by HERA data. Uncertainties due to different extrapolation 
approaches are about 30%, and a factor of two at E u = 10 12 GeV if the gluon saturation effect is 
taken into account [54] (see e.g. [12] for a brief discussion on this). 

Beyond the SM, the uncertainties in the neutrino interaction are conceivably larger. In many 
extensions of the SM, neutrinos can interact with nucleons via additional channels, the rates for 
which exceed the SM one largely [19]. If neutrinos become comparably strongly interacting as 
protons at an energy around 10 ia6 GeV, ultrahigh energy cosmic neutrinos could have already 
manifest themselves as the observed cosmic ray events beyond the Greisen-Zatsepin-Kuzmin 
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Figure 2: Predictions for the cosmogenic neutrino flux per flavour. Solid line: flux from [14] 
(FKRT) with the proton spectrum as the best fit to the Akeno + AGASA cosmic ray data, 
corresponding to a = 2.57, n = 3.30 and E max = 10 125 GeV. Dashed line: flux from [48] (PJ) 
assuming a maximum energy of E miLX = 10 12 ' 5 GeV for the ultrahigh energy cosmic protons. Long 
dashed line: flux of electron antineutrinos from neutron /3-decay (parent heavy nucleus 56 Fe, see 
Appendix). Also shown are the differential upper limits on the ultrahigh energy neutrino flux 
per flavour. Lower solid line with squares: upper limit by RICE, Eq. ()17j). Upper solid line with 
triangles: the combined Fly's Eye + AGASA upper limit obtained in Ref. [11]. The dashed-dotted 
line and the hatched area are the sensitivity (defined as one event per year and per energy decade) 
of PAO for v e and for Earth-skimming v T [36], respectively. 

(GZK) cutoff. The idea of using this "strongly interacting" neutrino scenario [13] to solve the 
GZK puzzle is supported by the observation that the cosmogenic neutrino flux well matches the 
observed ultrahigh energy cosmic ray spectrum above the GZK energy Eqik ~ 10 10 9 GeV (see 
e.g. FKRT in Ref. [19] for a statistical analysis of this scenario). 

The neutrino-nucleon inelastic cross section at ultrahigh energies is constrained by several con- 
siderations. As s goes to infinity, the unitarity bound [55] limits the cross section to grow at most 
as <r tot < C ■ In 2 s. However, without the knowledge of the constant C this bound cannot be of 
practical use. On the other hand, the neutrino-nucleon cross section at high energies is related to 
the low-energy neutrino-nucleon elastic amplitude through dispersion relations [56]. Laboratory 
neutrino-nucleon fixed-target scattering experiments at relatively low energies can therefore indi- 
rectly observe or constrain anomalous enhancements of <j v n at high energies. The non-observation 
of ultrahigh energy neutrino-induced events reported by several experiments such as Fly's Eye, 
AGASA and RICE imposes upper bounds on the neutrino-nucleon inelastic cross section as well. 
They depend, however, sensitively on the neutrino flux. 

In what follows we derive upper limits on the neutrino-nucleon inelastic cross section by exloiting 
the predictions for the cosmogenic neutrino flux discussed in Sec. Eland the formulae presented in 
Sec. 121 First we note that Eq. (j2J) and Eq. (j5J) can be used to constrain new physics models which 
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E v [GeV] 


FKRT [14] 


PJ [48] 


10 iu 


1.3 x 1CT 1 


2.4 x 1(T 2 


10 10.5 


no 


1.4 x 10" 1 


10 11 


no 


no 



Table 1: Upper bound on the neutrino- nucleon cross section (in [mb]) derived from the AG AS A 
Collaboration search results [7]. 



predict an enhancement of the SM neutrino-nucleon cross section above some energy threshold. 
The "strongly interacting neutrino" scenarios are also subject to this constraint. If the jump of 
the neutrino-nucleon cross section is not strictly a step function, neutrinos with energies in the 
intermediate range, where the cross section is ~ 0.1 fib - 0.5 mb, can always initiate showers deep 
in the atmosphere [57]. 

Now we derive model-independent upper bounds on a^(E u ) from AGASA's search result on 
quasi-horizontal air showers [6,7]. From Eq. (J2J), we demand [11,12] 

A p t (P(E Bh )} (E V F U (E U )) (att(£„)> < 3.5/A, (14) 

in a sufficiently small interval A, where a single power law 

P(E Bh ) F V {E) att(£) oc BP , (15) 

is valid. The choice A = 1, corresponding to one e-folding of energy, is reasonable. In Table Q we 
present our results by exploiting the cosmogenic neutrino fluxes discussed in the previous section. 
We have assumed that the neutrino-nucleon interactions do not distinguish between different 
flavours, and that the total neutrino energy goes into visible shower energy, i.e. E s ^ = E u . 

We found that these bounds are applicable only for cx*^ < 0.5 mb. If neutrinos interact more 
strongly, they would induce air showers high in the atmosphere, thus avoid the bounds derived 
from the non-observation of quasi-horizontal air showers by ground arrays. They would instead 
contribute to the vertical showers. Therefore, for energies E v > 10 11 GeV, the AGASA search re- 
sults on the horizontal showers cannot constrain the neutrino-nucleon cross section, if the neutrino 
fluxes are at the level as the cosmogenic ones we exploited. 

Next we derive upper bounds on the neutrino-nucleon cross section from RICE's search results 
[8,9] using Eq. (jSJ. We have checked again that, in the interval A = 1, the single power law 
approximation F V {E) a^(E) V(E) atten(£ l ) oc E a is valid, where we define 

r90° 

atten(E) = / d6 2 vr sin 9 e m p . (16) 

We demand 

t'pu 



o- 



(E U F V (E V )) (a^(E u )) (attend)) (V(E sh )) < 3.09/A . (17) 



We use the data taken in 1999, 2000 and 2001 [9], where the total lifetime is t' = 3500 hrs. The 
upper bounds on the neutrino-nucleon cross section from the RICE Collaboration search results 
are listed in Table El and given in Fig. EJ For comparison, one of the predictions for the SM total 
cross section [51] is also shown. We found that the bounds derived from RICE's search result are 
applicable for a l u % < 1 mb. 
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E v [GeV] 


FKRT [14] 


PJ [48] 


10 iu 


1.2 x 1(T 3 


2.8 x 1(T 4 


10 10.5 


3.6 x 1(T 3 


7.2 x 10" 4 


10 11 


1.9 x 1(T 2 


3.8 x 1(T 3 



Table 2: Upper bound on the neutrino- nucleon cross section (in [mb]) derived from the RICE 
Collaboration search results [8,9]. 
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Figure 3: Model-independent upper bounds on the neutrino- nucleon inelastic cross section de- 
rived from the RICE Collaboration search results [7] , by exploiting the cosmogenic neutrino flux 
estimates by FKRT (solid line) and PJ (dashed line joining solid squares). To give an idea of the 
scaling with the inelasticity parameter, the dashed line joining the open squares (PJ) indicates 
the upper bound for (y) = 0.5. The dashed-dotted lines indicate the sensitivity (95% CL, for 
al°N < 4 mb) of PAO in 10 yr of operation assuming zero events observed above SM background 
(circles PJ, triangles FKRT). For comparison, also shown is the SM total (charged current and 
neutral current) neutrino-nucleon inelastic cross section [51]. 

5 Sensitivity of PAO to Anomalous Neutrino Interactions 

In this section we estimate the potential of PAO to probe physics beyond the SM. Given the 
prospects for fairly high statistics, detailed analyses of high energy neutrino interactions are in 
principle possible. In particular, if no enhancement of deeply developing showers is observed, 
PAO will be able to set stringent bounds on the neutrino-nucleon inelastic cross section. 
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Event rates 


FKRT [14] 


PJ [48] 


+ V e + 


0.06-0.09 


0.16-0.25 


V T 


0.19-0.54 


0.56-1.59 



Table 3: Yearly neutrino event rates expected at the Southern site of PAO, for the cosmogenic 
neutrino flux estimates by PJ and by FKRT. Numbers in the first line are for hadronic showers 
falling within the array and all (hadronic) showers, respectively, induced by all neutrino (+ 
antineutrino) flavours in the range 10 9 GeV < E u < 10 11 GeV. Numbers in the second line are 
for Earth- skimming tau neutrinos in the range 10 9 GeV < E v < 10 11 GeV, obtained by assuming 
strong energy loss due to deep inelastic scattering (DIS), and no DIS-loss [36]. 



In order to estimate these bounds, we again adopt the cosmogenic neutrino fluxes shown in 
Fig. |21 and assume that only SM sources of deeply developing showers are observed. Note that 
in contrast to AGASA, SM neutrino interactions lead to observable rates at PAO. In Table El we 
list the yearly SM neutrino event rates to be expected at the PAO, for the cosmogenic neutrino 
fluxes we reviewed in Sec. El To obtain these event rates we used the SM cross section estimate 
given in Fig. 01 [51]. 

Before proceeding, we verify whether hadronic showers may be a significant background in our 
analysis. The number of cosmic ray showers expected to be detected by PAO in the angular bin 
6 E (75°, 90°) is given by 

m =a L cos9 dn j E p(e,q) e3 F;hs{E) ^ ' (18) 

where F° hs (E) is the incoming flux of cosmic rays and E% and E2 are the minimum and maximum 
energy under considerations. For E\ = 10 10 GeV, the PAO detection efficiency P(E, Q) reaches 
100%, and so Eq. (fTK|) can be rewritten as 

^^AAQ (E*F^(E)) . (19) 

Now replacing in Eq. (fT$|) the observed isotropic flux in this region, (E 3 F° hs (E)) = io 24 - 5±0 - 2 eV 2 
m~ 2 s _1 sr _1 [58], we obtain dN p /dt ~ 317 yr -1 . A detailed background event estimate requires a 
convolution with the detector resolution, and it is beyond the scope of this paper. However, the 
(A max ) distribution obtained [59] through Monte Carlo simulations indicates that the probability 
a proton-induced shower with 10 10 GeV < E < 10 11 GeV leads to A max > 2500 g/cm 2 is < 10~ 4 , 
hence hereafter we neglect the hadronic background in our analysis. 

With this in mind, as a very conservative estimate of background events, we consider the expected 
SM neutrino showers with 10 9 GeV < E v < 10 11 GeV, given in Table 01 In Table H] we list the 
sensitivity of PAO (95% CL corresponding to 3.54, 4.24 events for FKRT and PJ, respectively [31]) 
for anomalous neutrino cross sections after 5 yr of operation. In Fig. El the 10 yr prospects 
(N = 3.96, 5.08, for FKRT and PJ, respectively) to improve these bounds are shown. Note 
that in 12 yr of running, the SM background of deeply developing showers at the Southern site 
will correspond to 3.09 events. In such a case, another technique to bound the rise of the cross 
section can be used. Namely, by separately binning events which arrive at very small angles to 
the horizontal and comparing event rates of deeply developing showers and Earth-skimmers, the 
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E v [GeV] 


FKRT [14] 


PJ [48] 


10 io 


1.0 x 10~ 3 - 1.9 x 10~ 3 


3.1 x 1(T 4 - 5.6 x 1(T 4 


10 10.5 


4.2 x 10~ 3 - 8.7 x 1(T 3 


1.1 x 10~ 3 - 2.1 x 1(T 3 


10 11 


2.7 x 10~ 2 - 7.8 x 10~ 2 


6.2 x 1(T 3 - 1.5 x 1(T 2 



Table 4: Sensitivity (95% CL) of PAO for the neutrino-nucleon cross section (in [nib]) derived 
assuming non-observation of deeply developiong showers above SM background in 5 yr of opera- 
tion. 

neutrino-nucleon cross section can be inferred [26,60]. This is because the flux of up-going r's 
per unit surface area produced by Earth-skimming neutrinos is inversely proportional to cx*^, 
whereas the rate of deeply developing showers due to neutrino interactions in the atmosphere is 
proportional to oy^- Therefore, we conclude that the full observatory (Northern and Southern 
sites) in 6 yr of running will reach the sensitivity to probe cross sections at the level of SM 
predictions. 



6 Conclusions 

The search for ultrahigh energy cosmic neutrinos is entering an exciting era. Recently, a wide array 
of projects have been initiated to detect neutrinos by searching for low-altitude quasi-horizontal 
showers, cascades in the Antarctic ice-cap, and radio emission from neutrino-induced showers. 
Some of these experiments, like PAO and RICE, have been taken data, and others, like the km 3 
IceCube telescope, are under construction. 

Although such a high energy cosmic neutrino flux has not been observed to date, it has long 
been known that there should be a guaranteed cosmogenic flux resulting from the interactions of 
cosmic protons with the CMB. Therefore, the non-observation of deeply developing showers or 
cascades in ice reported by neutrino-detection-experiments in conjunction with these cosmogenic 
neutrinos can be used to constrain the behaviour of the neutrino-nucleon inelastic cross section 
in the energy regime far beyond the reach of terrestrial colliders. 

In this work we derived model-independent upper bounds on the neutrino-nucleon inelastic cross 
section from existing search results reported by the AGASA and the RICE collaborations. The 
bounds apply for neutrino-nucleon cross sections smaller than 0.5 mb (AGASA), or 1 mb (RICE). 
Interestingly, the search result reported by the RICE Collaboration improves the upper bounds 
derived from the non-observation of deeply developing showers at AGASA by more than one order 
of magnitude. Therefore, in the presence of additional neutrino-emitting-sources, there is only 
little room for new physics contributions to the inelastic cross section. 

We have also estimated the sensitivity of PAO to neutrino interactions. This hybrid detector 
will facilitate powerful air shower reconstruction methods and control of the systematics errors 
which have plagued cosmic ray experiments to date. Through a comparison of deeply developing 
showers and Earth-skimming events, in 6 yr of operation the observatory will be able to probe 
neutrino interactions at the level of SM predictions, providing a final verdict on the rise of the 
neutrino-nucleon inelastic cross section at ultrahigh energies. 
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Appendix 

If ultrahigh energy cosmic rays are heavy nuclei, the relic photons can excite the giant dipole reso- 
nance at nucleus energies E> 10 11 GeV, and thus there should be accompanying photo-dissociated 
free neutrons, themselves a source of /3-decay antineutrinos. The decay mean free path of a neutron 
is cT n T n = (E n /10 u GeV) Mpc, the lifetime being boosted from its rest frame value r n = 886 s 
to its lab value via T n = E n /m n . Compared to cosmic distances >100 Mpc, the decay of even the 
boosted neutron may be taken as nearly instantaneous, and thus all free neutrons are themselves a 
source of /3-decay cosmogenic antineutrinos. The neutron emissivity C n (E n ), defined as the mean 
number of neutrons emitted per co-moving volume per unit time per unit energy as measured at 
the source can be estimated as follows. Neutrons with energies above 10 9 3 GeV have parent iron 
nuclei with T > To ~ 2 x 10 9 which are almost completely disintegrated in distances of less than 
100 Mpc [61]. 2 Thus, it is reasonable to define a characteristic time r r given by the moment at 
which the number of nucleons is reduced to 1/e of its initial value A, and presume the nucleus, 
emitted at distance d from the Earth, is a travelling source that at D w (d — cr r ) disintegrates 
into A nucleons all at once [62]. Then, the number of neutrons with energy E n = Ea/A can 
be approximated by the product of 1/2 the number of nucleons generated per nucleus and the 
number of nuclei emitted, i.e., C n (E n ) = NCa, where N = A — Z is the mean neutron number 
of the source nucleus. Now, to obtain an estimate of the diffuse antineutrino flux one needs to 
integrate over the population of nucleus-emitting-sources out to the horizon [63] 

MEv) = = -i- fdE n C n (E n ) [l-ex P f-§^)l f Q } de v ^{e v ) 

dE v 4itH J L V E n T n J J J Q de v 

J i S [Eh - E n e F ( 1 + cos 9 V ) /mj , (20) 



x 



where the r 2 in the volume integral is compensated by the usual 1/r 2 fall-off of flux per source. 
Here, H is the Hubble constant, E v and E n are the antineutrino and neutron energies in the lab, 
9 V is the antineutrino angle with respect to the direction of the neutron momentum in the neutron 
rest-frame, and e v is the antineutrino energy in the neutron rest-frame. The last three variables 
are not observed by a laboratory neutrino-detector, and so are integrated over. The observable 
Ey is held fixed. The delta-function relates the neutrino energy in the lab to the three integration 
variables. The parameters appearing in Eq. ()20j) are the neutron mass and rest-frame lifetime 
(m n and T n ). Finally, dP/de-p is the normalised probability that the decaying neutron produces 

2 Because of the position of 56 Fe in the binding energy curve, it is generally considered to be a significant end 
product of stellar evolution, and indeed heavy mass nuclei are found to be be much rare in the cosmic radiation. 
Thus, here we adopt iron as the pivot nucleus species, i.e., A = 56. 
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a V with energy e v in the neutron rest-frame. Note that the maximum 17 energy in the neutron 
rest frame is very nearly Q = m n — m p — m e = 0.71 MeV. Integration of Eq. (J2Uj) can be easily 
accomplished, especially when two good approximations are applied [64]. The first approximation 
is to think of the /3-decay as a 1 — > 2 process of 5m^ — ► e~ + V, in which the neutrino is produced 
monoenergetically in the rest frame, with e v = €q ~ 5mjv(l — m 2 /5 2 mAr)/2 — 0.55 MeV, where 
Sttin ~ 1.30 MeV is the neutron-proton mass difference. In the lab, the ratio of the maximum V 
energy to the neutron energy is 2e /m n ~ 10 -3 , and so the boosted V has a spectrum with Ey E 
(0, 10 -3 E n ). The second approximation is to replace the neutron decay probability 1 — e~ Drrin l EnTn 
with a step function 0(E™ x -E n ) at some energy E™ x ~ 0(D m n /r n ) = (D/10 Mpc) x 10 12 GeV. 
Combining these two approximations we obtain 



~dE v 



dE A 



8 vr e Hq J E w Ea/A 



(21) 



where E™ m = max{£?^r 



2t 



}, and E™ ax is the energy cutoff at the nucleus-emitting-source 



< A(D/ 10 Mpc) x 10 12 GeV. For C A oc E~ A a \ integration of Eq. flZU leads to 



IK 



{E v ) w 10 6 



/rimax 



E- 



max \ a 



v 



E TJ 



dJ^A 



dE, 



(22) 



where £^>10 6 - 3 (56/ A) GeV and 



El 



2e £7; 



max 
A 



A 



10 



7.3 



56 
T 



77imax 

E A 
10 12 GeV 



GeV. 



(23) 



The sub-PeV antineutrino spectrum is flat as all the free neutrons have sufficient energy i? n >.Er /A 
to contribute equally to all the V energy bins below 10 6 GeV. Taking a = 2 as a reasonable ex- 
ample, and inputting the observational value E AFq dJ 7 A >s /dE A \ r ~ 10 5 eV m~ 2 s _1 sr _1 [58] 
Eq. (|22J becomes [63] 



2 d^jy_ 

v dE^ 



(E v ) w 4 x 10 1 



56 
T 



eVm 2 s 1 sr~ 



(24) 



Note that the /3-decay process gives initial antineutrino flavour ratios 1:0:0 and Earthly 
ratios nearly 3:1:1. Compared to full-blown Monte Carlo simulations [65], this paper-and- 
pencil calculation underestimates the flux by about 30%. Of course the situation described above 
represents the most extreme case, in which all cosmic rays at the end of the spectrum are heavy 
nuclei. A more realistic guess would be that the composition at the end of the spectrum is mixed. 
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